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Abstract — The popularity of opportunistic spectrum access 
(OSA), implemented by cognitive radios (CRs), necessitates that 
these networks should be able to provide service differentiation 
for different classes of traffic. One of the classic schemes for 
such a goal is through priority queueing. As link interruptions 
in OSA networks are frequent, we discuss in this paper four 
priority queueing disciplines in the presence of interruptions: 
preemptive-resume, non-preemptive, exceptional non-preemptive 
and preemptive in case of occurrence of an interruption. Ana- 
lytical results, in addition to simulation results to validate their 
accuracy, are provided. 
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I. Introduction 

The deployment of wireless services in diverse applications 
and networks implies that wireless links may carry different 
types of traffic with variable requirements and importance. 
Meanwhile, it is expected that opportunistic spectrum access 
(OSA) implemented by Cognitive Radio (CR), due to its com- 
prehensive and powerful features, will be widely deployed in 
next generation wireless networks |[T|, |[2|. Thus, as illustrated 
in Figure [T] a cognitive radio link may carry different classes 
of traffic where it is obliged to frequently release the channel 
to primary users. 

It is therefore crucial to develop a model which enables us 
to analyze the traffic metrics, such as throughput and delay, 
for differentiated services in CR networks where the operating 
channel is experiencing frequent failures due to the appearance 
of primary users or quality variations. 

One of the classic schemes to provide service differentiation 
in communication networks is through priority queueing, 
which is well discussed in the literature. However, in CR 
networks, the channel (i.e., the server of the queue) is subject 
to interruptions thus priority queueing should be studied in 
the context of a queueing model with server interruptions. 
In this paper, we introduce a priority queueing approach for 
a CR link with multiple classes of traffic, where the link is 
subject to recurrent failures and interruptions. The occurrence 
of interruptions represent the time where the operating channel 
becomes occupied by primary users or its quality is unaccept- 
able due to other reasons such as deep fading. The motivation 
for this work is that not only classic disciplines of preemptive- 
resume and non-preemptive, which are defined based on the 
interaction of different classes of traffic, should be revisited 
in the presence of interruptions, new priority disciplines can 
also be defined considering the occurrence of interruptions. 
These new disciplines take into account also the interaction 



Fig. 1. A CR link with multiple classes of traffic. 



of interruptions and different traffic classes. We discuss two 
such queueing disciplines in this paper; so totally, four priority 
disciplines will be investigated: 

• Non-preemptive (Non) where the server can not be pre- 
empted when any traffic is in service (even in the first 
arrival recovery period); 

• Exceptional Non-preemptive where the server can not be 
preempted in general, but it can be preempted in the 
arrival recovery period, before the start of the real service 
of low priority packets; 

• Preemptive (Pr) where a high-priority class can preempt 
the server any time; 

• Preemptive in case of failure where the server can not 
be preempted during the service, but if an interruption 
occurs, after the interruption, the high priority packet (if 
any) will be served first. This scheme, to the best of our 
knowledge, is a new scheme which is very well applicable 
to CR networks. 

To tackle these priority queueing disciplines, we first dis- 
cussed in |3| a general queueing model (with single class 
of CR traffic) for a CR link where the link was subject to 
interruptions. The results of that work will be extended here 
to priority queueing for multiple classes of CR traffic. Note 
that even the results for the CR queue with a single class 
of traffic, which was discussed in (3) for limited cases and 
are recalled here with more discussions, are more general 
compared to other papers who model the single-traffic CR 
queue with interruptions as a priority queue with primary 
and secondary users ||4|, |[5|. The reason is twofold: First, 
priority queueing models will be very complicated when the 
distribution of the duration of the primary users' busy periods 



are general. Second and more important is that, interruptions 
in our model may represent the intervals in which, the user is 
performing spectrum handover (multi-channel scenarios), not 
necessarily the busy periods of the primary users. Although 
previous models can provide accurate results for simple cases, 
they do not address the general interaction between primary 
and secondary users which can limit their usefulness. For 
instance, their results can not generally be used to discuss 
the relation between different CR traffic classes, which is the 
objective of this paper. 

In the literature, the notions of queueing with interruption 
and priority queueing have been very well discussed [6]- 
pO| . In these models, the interruption periods are the busy 
periods generated by higher priority of traffic and the priority 
discipline is preemptive. However in this paper, we consider 
a general distribution for the interruptions and also we dis- 
cuss other priority disciplines such as non-preemptive (Non), 
exceptional non-preemptive (ENo) and preemptive in case of 
failure (FP). 

The results of this paper provide utility relations for the 
purpose of differentiated packet-level performance evaluation 
in cognitive radio networks. The interruption periods can 
model the time intervals in which the CR network is searching 
for a new channel. In this case, these periods can be called 
recovery periods |TT| . We focus on this model thus recovery 
and interruption periods are used interchangeably. Any new 
proposed recovery scheme in physical or MAC layer can 
be easily compared with previous ones just by deriving the 
recovery time distribution and changing the parameters of the 
interruptions in the queueing model. 

Moreover, derived relations can be employed in optimiza- 
tion models to adjust the performance of the traffic. In CR 
networks which are equipped to learning and decision-making 
capabilities, this aspect could be very interesting. For instance, 
proposed priority schemes can be employed in an optimization 
and decision-making problem where the objective is optimiz- 
ing some performance metrics for different classes of traffic, 
where the decision to be made periodically is the priority 
queueing scheme to be selected among these four proposed 
disciplines. 

The main contributions of this paper can be listed as follows: 

• A comprehensive queueing model for opportunistic spec- 
trum access (OS A); 

• Analytical utility relations for four priority queueing 
disciplines with interruptions; 

• A new priority queueing discipline which is preemption 
in case of failure; 

The reminder of the paper is organized as follows. In Section 
[Il| an M/G/1/ queue with interruptions and with a single class 
of CR traffic is discussed. The results are then used in Section 
[In| to solve four priority queueing disciplines in the presence 
of interruptions. Simulation results are presented in Section IV 
and finally Section |V| concludes the paper with some remarks 
on future research directions. 

II. General Queue (Single Traffic Class) Analysis 

Consider a Cognitive Radio (CR) node (user) operating in 
a CR network using opportunistic spectrum access (OS A). 



This user continues operating over a channel, called its current 
operating channel, until it finds the channel unavailable due 
to the return of primary users or with an unacceptable trans- 
mission quality due to deep fading, for example. In a multi- 
channel scenario, the CR user starts a recovery process with 
a random duration to find a new channel and then restarts the 
operation over this new channel. For a single-channel scenario, 
the CR user waits for the channel to become available again. 
By modeling the event of missing or vacating the channel 
as a failure |11|, we can denote the operating (availability) 
and recovery periods as time to failure and time to recovery, 
respectively. 

Let two random variables Y and R represent the length of 
the operating and recovery periods with cumulative distribu- 
tion functions Fy and Fr, respectively. As illustrated in Figure 
|2] the operation of this user can be modeled with alternating 
renewal processes, assuming homogeneous channels. From the 
queueing point of view, the output buffer of the user can 
be modeled as a queue with random service interruptions. 
Random interruption implies that both the time of the failure 
occurrence which triggers a recovery, and the duration of 
the interruption (recovery time or waiting time) are random 
variables. If we assume that the arrival process of all classes 
of traffic follows a Poisson distribution with a general service 
time distribution, the operation of the node can be modeled as 
an M/G/1 queue with random service interruptions. 

As can be seen in Figure [2j we are now ready to mathemat- 
ically tackle the queueing model where the random variable 
Y represents the availability (ready for use) periods and R 
stands for the interruption (recovery) periods. Let us define 
also a new random variable C which is the sum of two random 
variables Y and R and represents the length of a cycle made 
of alternating renewal processes Y and R. The distribution of 
C can be given from the convolution of the distributions of Y 
and R. 

The service rate is fixed, so the distribution of the packet 
length defines the distribution of the real service time of the 
packets. The fixed service rate may imply however a fast 
varying channel whose average service rate is used as a fixed 
service rate in this model. The random variable T is used 
to represent the real service time (transmission time) of the 
packets. We call it real service time because the time that a 
packet spends in service in such a queue with interruptions 
includes both the real service time and the interruptions that 
may occur during its service. The notion of completion time 
is used for the whole time in service. A represents the Poisson 
arrival rate. The queue size is assumed infinite, so the main 
performance metric is the total time spent in the system (called 
system time or sojourn time) and in the queue (queueing 
time), represented by D and W respectively. Throughout the 
paper, for any random variable Z, fz{.) and Fz{.) respectively 
represent the probability density/mass function (PDF or PMF) 
and the cumulative distribution function (CDF) of the random 
variable Z. Moreover, Z{s) represents the Laplace- Stieltjes 
(LST) of the distribution of the random variable Z (Fz{.)). 
All the parameters and notations are listed in Table |Il| 

As can be seen in Figures [3] and |4] and discussed in J3J, 
we can distinguish three types of packet in such a queueing 
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Fig. 2. Queue Model with Interruption. 

model: 

1) When the packet enters an empty system and the server 
(i.e., current channel) is available, its real service starts 
immediately (See Figure [3] 1. The index 'a' is used to 
designate this case); 

2) When the packet enters an empty system during an 
interruption period such that the server is not available 
and a recovery before starting serving the packet needs 
to be completed (See Figure [3]2. The index 'u' is used 
to designate this case); 

3) When the packet enters a busy system and is queued, 
its service starts immediately after the completion time 
of the previous packet (See Figure |4] The index 'b' is 
used to designate this case). Note that in this case, the 
completion time of the packet is always started within 
an availability period (similar to the case 'a'). 

As there are two perspectives for queueing with interruptions, 
we discuss both in the beginning. Their difference is for the 
packets which arrive during an empty system and the system 
is in the recovery period (case 'u'). 

In the first model, the arrival time is the start of the 
completion time in the case 'u' (Figure |3] 2). However, the real 
service time of a packet is started when the current recovery 
period is completed. So, the remaining time of the current 
recovery period (Rr) is not accounted as waiting time, but the 
completion time. 

In the second model which is used more in the literature, 
the start of the completion time is always the start of the real 
service. Therefore, the remaining time of the recovery period 
in arrival for the case 'u' (Rr) is considered as waiting time. 
To distinguish these two, we use a star, '*', for the first case. 
Note that as expected, the system time for both cases will be 
the same and the difference is for the completion time and 
waiting time of the case 'u' packets. 

It is worth mentioning that when the transmission of a 
packet is interrupted, the whole packet should be generally 
re-transmitted and the assumption of a partial transmission 
(resume of transmission) is not realistic in some applications. 
This case will be postponed for future and we discuss here 
only the cases that partial transmission is possible and thus 
after an interruption, the node transmits the remaining part of 
the packet (service-resume). 

A. Case 'a': Arrival to an Empty -available System 

For packets of case 'a', the waiting time is zero. As the 
first availability period is started in a random point during 
an availability period, we assumed that the remaining time 
until the next interruption is distributed with a random variable 
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Fig. 3. Completion time for the case 'a' (1.) and case 'u' (2.). 

called Yia (first availability period for case 'a'). As illustrated 
in Figure [3] 1 , the completion time of the packet in the first 
case, Xa, can be given by: 



--R 



Yia>T 

+ R1+Y2 + R2 + ... + Yk + Rk + Tr, Otherwise. 

(1) 

where Tr stands for the real transmission time of the last part 
of the packet in the last operating period. For a service time 
distributed with the random variable T, we have: 



Yta 



Yk+Tr. 



(2) 



That is, the sum of transmission times in each availability 
period should be equal to T. The unknown parameter is 
K which is the number of availability periods during the 
completion time. If we consider only the operating periods, 
Yi, as a renewal process, K in all scenarios is the number of 
renewals of Y during the real transmission time of a packet 
where its distribution can be found from the renewal theory 
results (I2I, (ig. 

The exact distribution of Yia depends on the time that the 
user has switched to this channel, i.e., after the last recovery 
period (multi-channel scenarios), or the time that the system 
has become empty. Finding the length of this time interval 
can be complicated for general distributions. If we find the 
distribution of Yi^, the moments of the number of renewals, 
K, can be given by: 



C[ma{t)] 



Yia{s) 



s(l-Y{s)) 



(3) 



where ma{t) is the mean renewal function (average number 
of renewals during (0, t]) composed of instances of Yia and Y 
|[T2|, | p3| . One can find ma(t) from the equation above and 
thereby derive E[Ka]\{T = t) [j Unconditional mean number 
of renewals can then be given by: 



E[Ka] = / ma(t)fT(t)dt (4) 

Now let rn^(t) be the second moment of the number of 
renewals in time t. The LST of m^(t) can be found as: 



C-[ml{t)] 



Y\a{s){l + Y{s)) 



(l-Y{s))' 



(5) 



so that E[Kl]\{T = t) is found. 
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Fig. 4. Completion time for the third case when the packet enters a busy 
system and is queued (case 'b'). 



When the first two moments of Ka is found, using Eq. ([T]), 
we can give the first moment of Xa by: 

E[Xa]=E[T]^E[Ka]E[R]. (6) 

The second moment can similarly be found equal to: 

E[Xl] = E[T^] + 2E[TKa]E[R] 

^E[K,]{E[R^] - {E[R]f) + {E[R]fE[Kll 

(7) 

using the fact that the number of renewals, Ka, is independent 
of the recovery process, but not the service time. The second 
term in the calculation of £^[Xa] is a random sum whose 
moments can be found easily q] E\TKa] can be found by 
conditioning on the value of T, as follows: 



E[TKa] 



Jo 



E[TKa]\{T = t)fT{t)dt 

tma(t)fT(t)dt. 



(8) 



B. Case '/?'; Arrival to a busy system 

If we find the number of renewals for the case 'b' equal 
to i^5, based on the distribution of the first availability period 
which is here distributed with another distribution, Y15 (Figure 
|4]), we can use all previous results just by replacing Yia with 
yi5 and Ka with i^^. However in this case, the waiting time 
is not zero. 

C Case 'u\' Arrival to an empty -unavailable system 

For the case 'u', as we discussed, we have two models. 
In the first model, the completion time of the user is started 
within a recovery period. As illustrated in Figure [3]2, we have: 

X: = i?, + Fi + i?2 + ... + Ik + Rk+1 + Tr 
X, = Fi + i?2 + ... + Ik + Rk+1 + Tr, (9) 

and 



w: = 0, 

Wu = Rr- 



(10) 



The condition is similar to Eq. ([2]), replacing Yia with 
Yi. Similar to the discussion for the first operating period in 
the previous case, here the distribution of the first recovery 
time {Rr) (remaining time of R) can be different from other 
interruption intervals because the server knows how long it 
has spent so far for the recovery. 

^For instance, Var{KR) = E[K]Var{R) + {E[R])'^Var{K) 



The number of renewals is independent of the model and it 
can be found by replacing Yia in Eqs. ([3]) to ([5]) by F because 
here the first availability period is also distributed as Y. The 
moments of the completion time can be found similarly with 
a small difference that we should add the first recovery period 
(with distribution R^) to the completion time of the first model. 
In other words, we have: 

E[X:]=E[Rr]+E[X^], 
E[Xl*] = E[Rl] + E[Xl] + 2E[Rr]E[X^]. (11) 

As discussed, we can write: 



^u = ^a when ^Ka^KXia^Y) 
= Xb when (^K,^K,Y^,^Y)' 

The distribution of Rr can be written as: 
Rr = R-A\{R> A), 



(12) 



(13) 



where A is the inter-arrival time. That is, the time to arrival of 
the first packet is an inter-arrival time but conditioned on the 
fact that it should be less than R. The moments of Rr based 
on R and A are found in Appendix. 

D. Overall Performance Metrics 

The first moment of the general completion time can be 
given by solving the following equation: 

^[X*] = p*^[X5] + (l-p*) {PaeE[Xa] + (1 - Pae)^K]) , 

(14) 
where p* equal to XE[X*] stands for the probability of system 
being not empty and P^e represents the average probability 
that the server is available when the system is empty. From 
|6, Eq. (23)], Pae can be found equal to: 



Pn. = 1 



(l-Fy(A))(l-F^(A)) 

XE[Y]{l-FY{X)FR{X)y 



(15) 



where Fz{.), for an arbitrary distribution function Fz, can be 
defined as |6^ P.753]: 



Fz{0) 



f 
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e-'^dFz{u). 



(16) 



For = X and Z = Y or R, this function gives the probability 
that the length of the operating and recovery periods is less 
than an inter- arrival time. As expected, when A is small so 
that the inter-arrival times are long compared to Y and R, 
several operating and recovery periods can be observed during 
an inter-arrival period. Therefore, for small A, E[R] and E[Y], 
we can write: 



E[Y] 



(For small A). (17) 



" E[Y]^E[R] 

In this case, since Pae is a function of E[Y] and E[R] and 
Y and R are assumed stationary processes, we can assume 
Pae independent of Y. Also in general and for the regular 
distributions, P^e is a function of only the moments of Y 
and R. Therefore, we assume for simplicity that there is no 
correlation between Pae and Y and R. Note that Pae is a 



conditional probability, conditioned on the fact that the system 
is empty. 

We also use the notation of X* for the completion time of 
a packet which arrived to an empty system (cases 'a' and 'u' 
together). Its expected value can be given by: 



The second moment of the general completion time can be 
found equal to: 



E[X'^]=pE[Xt]^{l-p)E[Xn 



(19) 



Now to derive the waiting time, we use the same approach 
as an M/G/1 queue. When a packet arrives, it waits for the 
remaining completion time of the packet in service (if any) and 
then the completion time of all packets in the queue. The point 
is that for the packets which are in the queue, the completion 
time is always distributed with X5 (they are queued, so they 
have not arrived to an empty system). However for the packet 
which is in service, the general completion time should be 
used because no knowledge is available to know whether this 
packet has been of case 'a', 'b' or 'u'. Thus, we can write: 



£;[w*] 



\E[X 



2*1 



2{l-\E[X^]y 



(20) 



We can see that the average remaining completion come for 
a packet in service (similar to the remaining service time in 
an M/G/1 queue) can be written equal to: 



E[J, 



c\ 



^^ '2E[X*] 



A 



E[X^ 



(21) 



E. Second Model of Case 'u ' 

We can all the parameters above for the the second model of 
case 'u'. E[X] and £'[X^] can be found similarly by replacing 
XI with Xu. As discussed, Jc and E[T] = E[W\+ E[X\ = 
E[W*] + E[X% remain the same. Therefore, E[W] can be 
found as follows: 



E[W] = E[W*] + E[X*] - E[X]. 



(22) 



F. Busy Periods 

Using the same idea of calculating the distribution of the 
busy periods in an M/G/1 queue without interruption fSl, we 
can find the busy periods of our queue with interruption. We 
know that the first completion time in a busy period is an 
instance of X^. However, for other busy periods which are 
initiated during Xg (the idea of serving in an LCFS manner 
(8J), the busy period is started with an instance of X5 because 
the packets enter a non-empty system. Therefore, we can find 
the LST of the busy periods equal to: 



B{s)=Xe{s^\-\B^{s)), 



(23) 



where ^^(s) is the LST of the busy periods which are initiated 
during X^ thus with an instance of X5. Therefore, ^5(5) itself 
can be found from the following equation: 



B^{s)=X^{s^\-\B^{s)). 



(24) 



From the equation above, we can find the first and the second 
moments of the Bi)(t) and from there, the first and the second 
moments of the general busy periods, equal to: 



^[""l - l-\E[X,\ 



(25) 



E[X:]=PaeE[X^] + {l-Pae)E[X*^]. (18) and 



E[B-^] = \E[Bl]E[X,] + (1 + \E[Bb]fE[X^]. (26) 

The moments of the busy periods which start with an instance 
of Xh, E[Bh] and E[B'^], can be found from solving Eq. (24) 
equal to: 



E[Bi, 



E[Xb 



and E[Bl 



m! 



(27) 



1-\E[X,] ' "' {l-\E[X,]y 

Similarly, we can define the busy periods which are ini- 
tiated with any initial delay. Note that similar to remaining 
completion time, the busy periods are the same in both models 
because even in the second model, the arrival of the packet 
should be considered as the beginning of the busy period. 

G. Approximation for the Distribution of Yia and Yn) 

To find Yia, as discussed in |6|, we should consider two 
possibilities. When the last packet before our packet in con- 
sideration, which has just arrived, has left the system in the 
same availability period, and the case that the last packet has 
left sooner, so during this availability period, the system has 
been empty. In the former case, we can see that the remaining 
time of the availability period after the departure of the last 
packet before us is distributed as Yi^. The time between its 
departure and the next arrival (of our packet in consideration) 
is Exponential with parameter A. Therefore, the remaining time 
of the availability period can be written as: 



Yia I (Last departure in the same availability period) 



(28) 



(29) 



= ri5-A|(yi5>A). 

In the latter case, we can simply replace Y15 with Y. We have: 
YiaKLast departure in previous availability periods) 

= Y -A\{Y> A). 

Finding the probability of occurrence of each case is cum- 
bersome and, as discussed in |[6|, needs some simplification 
assumptions for Yi^. As the distribution of Yia and Yi^ is 
complicated, we can follow two approaches. The first one is 
finding lower and upper bounds, as discussed in |[6| Eq. (5)], 
for yi5 and consequently for the general completion time. 

The next approach is proposing an approximation for Yia 
by assuming that Yia may be started uniformly during an op- 
erating period which is sometimes called random modification 
of Y \1\ or equilibrium excess distribution |6|. We can write 
the distribution of Yia based on Y equal to: 

l-Fy(t) 



fnAt) = 



E[Y] 



(30) 



The n-th moments of Yia based on the moments of Y can be 
written as [12J : 



E[yra] = 



{n^l)E[Y]' 



(31) 



Note that this is an approximation for Yia', thus, it is expected 
that we observe some difference between the simulation and 
analytical results. By this assumption, we approximate that 
X5 has almost the same distribution as Xa. That is, ^^[Xfo] ^ 
E[Xa] and E[X^] ^ E[X^]. Note that by this assumption 
for Yi5, the probability of the occurrence of each one of the 



possibilities for Yia (Eqs. ( [28] ) and ( |29| )) can be found equal 
to|6j: 

Pr (Arrival in the same availability period) 

(32) 



the recovery period in the arrival point of the packet which 
initiates a busy period is modeled as vacation of the server 
to use the available results in the literature for M/G/1 queues 
with vacations (queues with initial setup time) |[8|. The initial 
setup time for a packet which initiates the busy period is Rr 
with the probability (1 — Pre) and zero otherwise. Thus, we 
can find the moments of the initial setup time based on the 
moments of Y and R. That is, E[S] = (1 - Pae)E[Rr] and 
E[S'^] = (1 - Pae)E[Rl]. From |8, Eq. 2.44a], we have: 



= 1 



Jo 



E[W] 



XE[X^ 



This approximation helps us to have an equilibrium renewal 
process for cases 'a' and 'b' fVT\ which is much easier to 
calculate. The LST of the moments of the number of renewals 
with this assumption can be written as: 

1 ... t 



2(1-A^[X5]) 



and 



2E[S] + XEIS"^] 
2(1 + \E[S]) 
E[R^] 
2{l-\E[Xr,]) ' 2{E[Y]+E[R]y 

E[D]=E[Xb]+E[W]. 



(39) 



(40) 



C[ma{t)] = 



and 



-E'[F]S'^ E\Y\ The steady-state probability of system being empty, Pq, can 

he siven hv: 



t^iml^t)] = 



Y{S) 



be given by 



(34) 



S'^E[Y]{l-Y{S))' 

The unconditional moments can then be found similar to Eqs. 
(|6J( and ^, but with the new distributions for ma{t)^mh{t) = 
=^|(T = t)andform2(t)«m2(i). 

H. Exponential Availability Periods 

As in the paper it is mostly assumed that availability periods 
are distributed with an Exponential distribution, we discuss 
this case separately here. Y is Exponentially distributed with 



Pn = 



^ 1 - XE[X,] 

E[I]^E[Bs] 1^XE[S] ' 



E[I] 



(41) 



the parameter a (i.e., Fy = Fy^^ = 1 



*) thus we have 



Yia = Yii) = Y and Ka = Ku = i^6 because the remaining 
time in an availability period is still the same exponential 
distribution. 

Using Eq. ([s]), m{t) can be found equal to at. From Eq. 
(pb, the second moment of the number of renewals, im?{t), is 
equal to rin?{t) = a^t^ ^at. These results verify our relations 
since the number of renewals for exponentially distributed 
availability periods is a Poisson process with the average at. 
Therefore, we have: 

E[Xa] = E[Xb] = E[X^] 

= E[T] + aE[T]E[R] = E[T]{1 + aE[R]), (35) 

/>oo 

E[TKa] = / at^Mt) = aE[T% (36) 

Jo 



E[X^] = E[T^]{1 + aE[T]f + aE[T]E[R% (37) 



^[Xf ] = E[X'J + E[Ri] + 2E[Xa]E[Rr]. (38) 

An interesting point when the availability periods are Ex- 
ponentially distributed is the fact that considering the second 
model assumed for case 'u', the completion time of all packets 
is the same (equal to X5). So, to find the waiting time, 
we can model our queue as a queue with vacation. That is. 



where E[I] = j is the average of Idle periods (no packet in 
the system), and E[Bs] is the average of busy periods initiated 
by S -\- Xi) which can be found from Section |II-F| 

III. Priority Queueing 

The results found so far can be used to tackle priority 
queueing disciplines over this general queueing model with 
Np classes of traffic. We define and apply four different 
preemption schemes: classic schemes of non-preemptive and 
preemptive -resume p4| , an exceptional non-preemptive and 
a new one which is called preemption in case of failure. A 
non-preemptive scheme permits a low priority (LP) packet 
to be totally transmitted before serving a new arrived high 
priority (HP) packet. This implies that even if an LP packet 
arrives in an empty system during the recovery, incoming HP 
packet in the same recovery period should also wait for the 
termination of the LP service. While in the exceptional non- 
preemptive scheme, the HP packet which is arrived in the same 
arrival recovery period starts the service. But, If LP starts the 
real service, the service can not be preempted. A preemptive 
discipline preempts the service from the LP packets in any 
time as soon as an HP packet arrives. The new scheme which, 
to the best of our knowledge, is a novel scheme and very well 
applicable to cognitive radio networks, does not preempt the 
service during the real transmission of the LP packet. But, if a 
failure occurs which necessitates a recovery, after the recovery, 
the HP packets (if any) will be served first. When the HP 
queue is empty again, the transmission of the LP packet from 
the same state is resumed. 

For simplicity of presentation, we will use parameters with 
an index 1 for high priority (HP) packets and with an index 2 
for low priority (LP) ones. We focus on two classes however 
it can be, sometimes, extended to multiple classes of traffic. 
Furthermore, in the rest of the paper, when notations A, p 



and A are used without any index, they represent combined 
parameters. That is: 

Np 

A = ^A„ (42) 

Np Np 

P = JlpJ=Jl>^J^[^i]^ (43) 

and 

A = min{Ai, ..., A^J ^ A - EXP{\). (44) 

We also use the notation of PaeW and i?^(A) to highUght 
that Pae and i?r in Eqs. ( 15 ) and ([T3|, respectively, should be 
calculated with combined A and A. 



A. Non-preemptive Priority 

In a non-preemptive scheme, when any packet from any 
class goes to service, its service will not be interrupted by other 
packets. This implies that the completion time of any packet 
of any class will be the same as the original M/G/1 queue with 
interruptions, as already discussed. However, even though the 
completion time of all cases is the same as the original queue, 
the average completion time for the first model (with '*') is 
different because the probability of occurrence of each case 
('a', 'b' and 'u') is different. That is, in a queue with multiple 
classes of traffic and when the service can not be preempted 
from the packet in service, the probability of arrival in an 
empty system is less. Thus, X5 plays a more determining role 
and consequently the average completion time in such a queue 
is less, compared to a queue with a single classes of traffic. 
However, the waiting time naturally increases. 

1) First Approach: If we return to the original queue and 
replace the probability of system being empty from \E[X*] 
to T!^Zi\jE[X*] in Eq. (14) for all classes of traffic, we 
will have Np equations and Np unknowns which enables us 
to find the new average completion time of the packets in this 

is replaced with ^ \j 



system. To find P^e, A in Eq. 



15 



for 



all classes of traffic. Consequently, the second moment can be 
found. Then, similar to an M/G/1 queue 1 14] we have: 

E[r] 



E[W-] 



(45) 



We used the fact that Cyr{T) = E[Xi)] for all queued packets. 
J* stands for the remaining completion time of the packet in 
service, which is the same for all priority classes. It can be 
given by: 



N^ 



A. 



^[^*] = Et^[^^-*'] 



(46) 



where Np is the number of priority classes. This implies that 
for the packet in service, no knowledge is available whether 
it has been a packet of case 'a', 'b' or 'u', so a the general 
completion time is used. However, the denominator represents 
the completion time of the queued packets which is X^ for 
class j. 

When Y and arrival processes are all exponentially dis- 
tributed, which is of our interest in this paper, we know 



that Xb = Xa = Xu. Thus, using the idea of Eq. ( [39] ), a 
closed-form relation can be given for the waiting time of the 
packets. For this aim, the queue can be modeled as a queue 
with an exceptional completion time for the first packet which 
initiates a busy period. As saw for the original queue, this 
exceptional service is represented by X* (Eq. ([19])). However 
here, X* should be updated by the combined probability of 
arrival during Y or R. We define: 

Xf = PaeWiXl) + (1 - Paeim^l + Rr{A)), (47) 



N^ 



^e = E^xr. 



i=i 



(48) 



First two moments of Xg can be found easily. Then, we can 
write I8J: 

(1 - p,){E[Xl^] + ^[X^]) + XE[Xl^]E[Xl] 



E[T^] = 



1 + A^[Xe] - p 
A[(l - p)E[X^,] + ^[Xe](Ef^i A,^[(X^')^ 



2(1 + A^[Xe] - />)(! - EU />.)(! - E;=i Pj) 

(49) 

The first term represents the average completion time and the 
second term, the average waiting time. Note that the equation 
above is a closed-form relation of the general case (provided 
in Eq. ([45])), for exponential availability periods. 



Similar to Eq. (41), the steady-state probability of system 
being empty, Pq, in this queue can be given by: 



Po 



1-p 



(50) 



-p^XE[Xe]' 

2) Second Approach: Another idea is to use available 
results in the literature for the mixed priority queueing dis- 
ciplines. As interruptions have a preemptive behavior, we can 
model the interruptions as the highest priority type of traffic 
whose inter-arrival times are distributed with a random vari- 
able Y and whose busy periods are distributed with a random 
variable R. From these two, we can have an estimate for the 
service time of these virtual highest priority packets whose 
activity models the interruptions. Closed-form relations can 
not be derived in general since for any distribution of R and 
Y, a different formula for the busy periods and consequently 
for the service time exists. If Y is Exponentially distributed 
with parameter a, we can assume an original M/G/1 queue 
for the interruptions. Therefore, from the following equation, 
which is the distribution of the busy periods of a regular M/G/1 
queue fS], we can find the first two moments of the service 
time of the virtual packets which form the interruptions: 



E[T,] 



E[R] ^[j..j 



l^aE[R] 



a-aR{s)), (51) 

= E[R^]{l-aE[T,]f. (52) 



where Ty stands for the service time of the virtual packets who 
form the interruptions. Then, the queue can be modeled as a 
priority queue where the highest class of traffic (virtual) be- 
haves preemptively, but other classes behave non-preemptively. 
From ( 1 8 1), the extensions of the P-K formula can be used 
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Fig. 5. The difference between Non-preemptive and Exceptional Non- 
Preemptive schemes. 



to find the moments of the waiting time of other classes of 
traffic. For the class i, it can be given by: 



E[W^] 



E[J] 



(1 - aE[n] - ... - A.-i^[T._i])(l - aE[n] 



where E[J] is equal to: 



■X^E[T^-l])' 

(53) 



1 "^^ 1 

E[J] = -aE[T^]^^-X,{E[Tf]). (54) 



j=i 



The system time can be written as: 

Em 



E[Di_ 



■ aEm 



+ E[W^] 



(55) 



Note that this approach provides an approximation because 
we know that the real distribution of the busy periods in an 
M/G/1 queue is a complicated function built on the Bessel 
function |8|. For other distributions of availability periods, if 
we are able to find the probability generating function (PGF) 
of the number of arrivals, we can apply the same approach 
to find the moments of the service time of the virtual packets 
which form the interruptions. 



B. Exceptional Non-Preemptive 

If we remember two different models which discussed in the 
original queue (with and without '*'), we can discuss a new 
non-preemptive scheme with a small difference that if before 
the start of the real service time of an LP packet, an HP packet 
arrives, the HP packet is served first. As illustrated in Figure 
[5] this case happens only when an LP packet arrives into an 
empty system and during the recovery period (case 'u' or the 
original queue), and an HP packet also arrives in the same 
recovery period. In other words, this the a non-preemptive 
scheme when the completion time is defined without '*' (start 
of the real service). 

Similar to the previous discipline, this queue can be modeled 
as a queue with initial setup time with the difference that 
during the setup time, the service is preemptive |8|. As the 
arrival process is Poisson for both classes of traffic, the setup 
time is the same for both classes and similar to the previous 



discipline. The service time of different classes of packet can 
be given by: 



>Ar„ 



E[W,] = 



E-ZrXjE[{Xi)'] 



{l-p){XE[S']^2E[S]) 

(2(1 + A^[^] - p))(i - e;=i p.)(i - e;=1 pj) 



E[D,]=E[W,]^E[Xl]. 

The probability of system being empty is equal to: 

1-p 



Pn = 



1 + XE[S] ' 



(56) 

(57) 

(58) 



C Preemptive -resume Priority 

In this scheme, the highest class is not affected by other 
classes of traffic. So, the completion time and the waiting time 
of the highest priority packets can be found from the original 
M/G/1 queue with interruptions. Let us find the completion 
time and the waiting time of other classes. 

1) First Approach: An approach is to use available results 
in the literature for the preemptive-resume schemes with 
multiple classes of traffic, similar to the one that we used 
for the previous scheme. As interruptions have a preemptive 
behavior and we are discussing the preemptive schemes, we 
can model the queue as a preemptive-resume queue with 
Np-\-l classes of traffic. Then the following famous extensions 
of the P-K formula for preemptive-resume schemes ( jSj) can 
be used to find the moments of the waiting time of other 
classes of traffic: 



E[W^] 



E[J^] 



(1 - aE[n] - ... - A._i^[T,_i])(l - aE[n] - ... - A.^[T,_i]) ' 

(59) 



where E[Ji] can be given by: 



E[Ji] = laE[T^] + J2l^AEm, 



(60) 



and the moments of Ty can be found from Eq. (52). 

2) Second Approach: Another approach, which is provided 
here because it will also be used for the third priority scheme, 
is to find the distribution of availability and unavailability 
periods from the perspective of the LP packets. Let us call 
them respectively Y2 and R2. Here, unavailability implies 
both due to the activity of HP users and due to interruptions. 
As illustrated in Figure [6]a, it is clear that Y2 is the minimum 
of two random variables Ai and Y which means between the 
interruptions and the HP packets, any one which arrives sooner 
initiates an unavailability period for LP packets. Ai stands for 
the inter- arrival time of the HP packets. Therefore, 



Y2 = min(r, Ai) ^ 1 - FY2{t) = (1 



■Fy{t)){l-FAAt))- 

(61) 
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Fig. 6. Availability and unavailability periods (Y2 and R2) from the 
perspective of LP packets. 



When Y is Exponentially distributed with parameter a, the 
distribution of Y2 can be given by: 



FY2{t) = l-e- 



-(Ai+a)t 



(62) 



To calculate R2, we should take into account which one of 
two events discussed above has occurred sooner. As can be 
seen in Figure [6] a, if the HP packet arrives sooner, the length 
of R2 is equal to one busy period of HP packets which is 
distributed with ^5. However, if an interruption arrives sooner, 
again two different cases may happen which discussed in 
sections b and c of Figure [6] During R, no HP packet arrives, 
so the length of R2 is equal to R (conditional), or at least 
an HP packet arrives. In the latter case, the length of the 
unavailability period from the perspective of the LP packets 
will be the interruption period which is extended with an HP 
busy period. We need the remaining time of the interruption 
period after the arrival of the HP packet. Recalled from Eq. 
(Tsl, this time is distributed sls Rr = R - Ai\{R > Ai) . After 



this time, an HP busy period is started with a duration equal 
to Bji^. Bji^ can be found from Eq. ( [23] ), replacing X5 by 
Xb^Rr. We have: 



R2= < 



'B, Pr{Ai<Y), 

(r\{R < Ai) PriY < AikR < A^)), 

\ai\{R > Ai) + Br PriY < A^kR > Ai). 

(63) 
From now on and for simplicity, we use the notation Zl^^^ 
for Z1\{Z1 < Z2) for any two random variables Zl and Z2. 
The probability of an HP arrival during R (i.e., Ai < R) can 
be given by: 



PaUR\{R = r) = l-e- 



-Air 



PaU 



(1 



-Air 



)dFR{r) 



For Exponential Y, we can write: 
Pr{Y < Ai) = 



a + Ai ' 



(64) 



(65) 



E[R2] 



Ai 



E[B, 



(1 - PaliR)- 



a + \i a + Ai 

-^PaUR{E[A<'']^E[BR^])]. 



-d/dXiR{Xi) 
R{Xi) 
(66) 



Using Eq.[92]for E[Rr] and Eq. 25 for busy periods, we have: 

E[Xl] ^ ^I^L _ J_ 



E[B 



rA 



E[Xl]^E[Rr] 
l-\iE[Xl] 



l-RjXi) 

XiE[Xl] 



(67) 

The second moment of R2 can be found similarly, as the 
second moment of the busy periods can be found from Eq. 
([26]) and the second moment of the conditional Ai and R can 
be derived from Eq. ( [9T] ) in Appendix. It should be taken into 
account that Br^ and Af^ are correlated, so E[BR^Af^] 
should be calculated separately, using, for instance, the same 
approach as in Eq. ([5]). 

From the equations above, one can find the moments of 
the availability and unavailability periods from the perspective 
of LP packets (Y2 and R2). Then, we return to the original 
M/G/1 queue with interruptions and replace Y and R in Eq. 
( [39] ) with Y2 and R2 respectively, to find the performance 
metrics of the LP packets. 

D. Preemption in Case of Failure 

In this priority queueing model, the type-1 packets can not 
preempt the service from a low priority packet when the latter 
is in transmission. However, the service can be preempted 
after an interruption. As expected, for both classes of traffic, 
the performance metrics are between two previous priority 
schemes. That is, for instance for the HP class, the average 
system time is more than the preemptive scheme but less than 
the non-preemptive scheme. Moreover, the completion time of 
the packets with the highest priority is not affected, so it can 
be still found from the original queue. However, the waiting 
time is affected. Other parameters which are the completion 
time of the LP packets and the waiting time of both HP and 
LP packets are discussed. 

1) Completion Time of the LP Packets: We want to find the 
completion time of the LP packets. The same idea used for the 
preemptive scheme can be employed with the difference that 
here, the HP packet can not preempt the channel before the end 
of the cycle. However, the difficulty is that the length of the 
new interruption periods from the perspective of LP packets, 
R2 ,depends on the remaining service time of the LP packets 
which naturally is not the same during the completion time 
of an LP packet. Unless, we assume a memory-less service 
time which is a special case and will be discussed later. The 
problem can not be therefore modeled as a renewal process 
because the instances of R2 are not identical. Thus, we provide 
some approximations. The approximations are discussed for 
two cases: when the availability periods are much larger than 
the service time of type-2 packets (Y » T2), called large 
scenarios, and when it is smaller (Y < T2), called small 
scenarios. 

Using the first assumption, it can be assumed that the 
service of a type-2 packet is finished in maximum two type-2 
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operating periods. The selection of two is a trade-off between 
accuracy and complexity. That is, if it is assumed that the 
service can be finished in only one operating period, the bound 
will be very loose, and if three type-2 operating periods are 
considered, the complexity increases. 

Assuming two type-2 operating periods is equivalent to ex- 
periencing at most one R2 interruption during the completion 
time. The duration of the interruption depends on the arrival 
of an HP packet, as already discussed. We can thus write the 
completion time of type-2 packets as: 




E[R2] 



Y>T2 
Y <T2. 



(68) 



A<^ 



where R2 stands for the length of the interruption period from 
the perspective of type-2 users. It can be given by: 

'R<^^ No HP arrival, 

R2 = <Bc^- {Y<^^ - Ai\{Ai < Y<^^)) HP arr. in Y<^^, 

Br HP arr. in R. 

(69) 

where Cr stands for the remaining time of the cycle after the 
arrival of an HP packet, and Bc^ represents the busy period 
which is initiated with Cr + -^6- However, if we return to 
Figure [6]a, we can see that the remaining time of F^^^^ which 
is represented by F<^2 —Ai, should be excluded from R2. In 
the third case, the HP arrival occurs in R. So, the busy period 
of HP packets is started with Rr + X5 and the length of the 
total interruption is Ai in addition to HP busy period. Thus, 
E[R2] can be given by: 



E[R2] = (1 - P,,,,< 



+ (Palic) 



Pae 



-d/d\iR{\i) 

R{\i) 

^E[R] + E[Y<^--Ai\[< Y<^^)]+E[Xl] 



(l-Ai^Kl) 
-E[Y^^^ - A^\{A^ <y<^^)]) 

E[Xl] + E[Rr 



+ {l-Pae){E[A^''] + 



{l-\iE[Xl]) 



(70) 



y<T2. 



where Paiic< is the probability of an arrival in C^ 
R, and Pae is calculated for HP packets. 

The second moment of R2 can be found similarly using 
the second moment of the busy periods and the relations 
provided in Appendix. However, again the correlation of 
random variables Br^ and ^4^ ^ should be taken into account. 

For the other approximation (small scenarios), as Y is 
supposed to be short, we assume that the duration of the HP 
busy periods is independent of the activity of LP packets. 
Therefore, the interruption periods from the perspective of 
LP packets can be assumed to have a similar distribution and 
thus, a renewal process can be considered. This implies that 
Eq. ([69]) is still valid however we eliminate the condition that 



Y <T. When a type-2 packet starts its service, as illustrated 
in Figure [7] it can hold the channel (available or unavailable) 
for a cycle and if there is an HP arrival during the cycle, it 
releases the channel to HP packets at the end of the cycle. So, 
the probability of releasing the channel to HP traffic in the next 
cycle is the probability of at least one HP arrival with the rate 
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Fig. 7. Cycles and holding periods for the LP packets in the discipline of 
preemptive in case of failure (FP). The LP packet can hold the channel for 
three cycles; however in a preemptive scheme (Pr), it can keep the channel 
for two cycles and releases the channel in the middle of the third cycle (HP 
arrival). 



Ai during the current cycle. We can find this probability equal 
to: 

Palic = Pr(HP arrival in C) 

{l-e-^nfc{c)dc, (71) 
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The interruption periods can be discussed similar to the 
preemptive scenario, as it is again assumed a renewal process. 
We have: 



R2 



R\{C<A,) 



No HP arrival in C, 



Ai I (C> Ai ) + Be, HP arrival in C. 



(72) 



where Cr stands for the remaining time of the cycle after the 
arrival of an HP packet. Due to the independence of Y and 
R, R\{C < Ai) = R\{R < Ai), which was found in Section 



III-C2| Similar to previous assumption, Cr itself is discussed 
separately for two cases: when the arrival has occurred in Y 
and when it has occurred in R, during the unavailability period. 
When the HP packet arrives in Y, the length of the re- 
maining time of the cycle can be given by the original cycle 
length, Cr = C = Y -\-R. Because, the remaining time of the 
availability period is the same as Y (memoryless property) and 
there is no condition for R. However, when the arrival occurs 
in R, we return to part (c) of Figure [6] That is, Cr = Rr- It 
can thus be written: 



d/dXiR{Xi) 



E[R2] = (1 - PaUc) 

+ {Palic) Pae{ 
-^{1 - Pae){E[A< 



R{Xl) 

E[R] 



■E[Y] 



(i-Ai^Wi; 



-E[Y]) 

■E[Rr], 



(1-Ai^K]) 



(73) 



The second moment can be found similarly. The moments of 
R2 then can replace R in Eqs. [35] and [37] to find the moments 
of the completion time of LP packets (£^[^2] and ^^[Xl]). 

2) HP Waiting Time: We can see that in this scheme, only 
HP packets of case 'a' (HP arrival during Y to the system 
which is empty of HP packets) are affected. For these packets, 
whose waiting time was zero in the previous case, the new 
waiting time will be a function of the LP queue. If the LP 



queue is empty, the waiting time is naturally zero. Otherwise, 
the waiting time is the minimum of the remaining time of the 
arrival cycle and the remaining time of the service of the LP 
packet in service. 

The difficulty is the dependency of both types of traffic on 
each other. That is, while the waiting time of the HP packets 
are affected by the lower class, both the waiting time and the 
completion time of LP packets are affected by HP traffic. This 
obliges us to suffice to approximations and bounds. 

For the exponentially distributed availability periods, which 
is the assumption in this part, we can see that the upper bound 
for the waiting time of HP packets is an M/G/1 queue with 
vacation (similar to the original queue) in Eq. (39)when S 
(initial setup time) is updated, as follows: 



(74) 



where Pae(Ai) and i^r(Ai) take only HP packets into account 
(not the total arrival rate), and Pq is the probability of system 
being empty of any type of packet. We can write: 
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XE[X^,] 
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\E[Xl 



\E[X^]) 2{E[Y] 
< E[Wi] < 

2E[S]^\iE[S^ 



(75) 



2(1 - \E\Xi, 



2{l ^\E[S]) 



This provides us an upper bound for the waiting time of the HP 
packets. As we will see later, by this bound we are assuming 
that HP packets should always wait for an LP packet (in case 
of arrival in Y) and the remaining LP service time is T2. 
Note that the probability of system being empty of any packet, 
Po, is the same for all priority queueing models discussed in 
this paper, regardless of the discipline. We have found this 



probability in Eq. (50). 

3) LP Waiting Time: When the moments of R2 and con- 



sequently X2 are known, a relation similar to Eq. (39) can 
be used to find the waiting time of LP packets (M/G/1 with 
vacation). To calculate the moments of vacation in this case, 
consider the arrival of an LP packet to a system which is 
empty of LP packets. If the arrival occurs in the time that an 
HP packet exists in the system, the vacation is the general busy 
period of HP packets, which is not the same as the original 
queue with single class of traffic, and difficult to derive. In 
case of arrival to a system which is empty of HP packets, 
two different cases should be discussed: An LP arrival in Y 
implies that the waiting time is zero. If it occurs in R, two 
cases may happen. No HP arrival occurs in the remaining part 
of R (called Rr), so the waiting time is Rr, and an HP arrival 
occurs, so the LP packet waits for the end of the HP busy 
periods. However, finding the moments of all cases discussed 
is difficult due to the dependency of HP and LP packets, and is 
not in the scope of this paper. We thus propose some bounds, 
as follows. 



Considering the Conservation Law (CL) in a queue with 
multiple classes of traffic fT5l, we know that 



n = XiE[Ti]E[Wi] + X2E[T2]E[W2] 



(76) 



is constant and the same for all priority disciplines proposed 
in this paper, regardless of the queueing discipline. In other 
words, as we found the waiting time of LP and HP packets 
in previous priority disciplines, we have the value of the sum 
above. If we find the waiting time of one of priority classes 
for the discipline of preemption in case of failure, we can find 
easily the other one. 

As an approximation for the first two moments of the 
completion time of LP packets has been found, the minimum 
waiting time of LP packets can be found by P-K relation, 
without considering the impact of vacations for the packets 
which enter an LP empty system. LP packets experience 
maximum waiting time in the preemptive-resume model. Thus, 
we can write: 



X2E[Xl] 



E[R^ 



E[R]) 



2(1-A2^[X2]) 2{E[Y] 
< E[W2] < 
E[W2] of preemptive scenario. 



(77) 



The Conservation Law helps us also to find another upper 
bound for the waiting time of HP packets, discussed in Eq. 
( [75] ). The point is that the maximum waiting of HP packets 
is when LP packets experience the minimum waiting time. 
Thus, by replacing the lower bound of ^[1^2] in Eq. (77) 
into Eq. (76) for E[W2], we can find another upper bound for 
HP packets. 



E. Exponentially Distributed Service Times 

The interesting point about the exponentially distributed 
service time is the fact that the remaining service time of the 
packets in each instant of time is still Exponential. This helps 
us to find more accurate relations for the last proposed scheme 
in which, only approximations were proposed. Discussion in 
1 14] indicate that the performance results with an Exponential 
service time can be considered as the worst case among 
scenarios with other common distributions for service time 
(e.g., Erlang or Deterministic). 

As we saw, the queue model for the HP traffic is an M/G/1 
queue with vacation where we found an upper bound for it 



in Eq. ( 75 ) and by using the Conservation Law. Exponential 
service time of LP packets enables us to find the moments of 
the vacation from the perspective of HP traffic. The vacation 
may also exist when the HP packet arrives in Y and should 
wait for end of the service of LP packets or an interruption. 



The HP waiting time can thus be given by Eq. (39) when S 
(initial setup time) is updated, as follows: 



S= < 



Rr{Xi) 



rp<Y 

-'-2 





(l-Pae(Ai)), 

Pr{Y <T2){PL\NH)Pa. 

Pr{Y>T2){PL\NH)Pa. 

Otherwise. 



(78) 
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where Pae(Ai) and i^r(Ai) take only HP packets into account 
(not the total arrival rate). E[S] can be written equal to: 



E[S] = (1 - Pae)E[Rr] 
+ {PL\NH)Pae 



a 



7 + a a + 7 



E[R]) 



■ (79) ^ 



The unknown here is Pl\nh^ which is the probability that 
the system is empty of HP packets, but there are LP packets 
in the system, in the point of arrival of an HP packet. The 
probability of system being empty of HP packets, Pq can be 
found from the original queue, replacing E[S] with the one 



that was calculated in (Eq. 79). Using these two, we can write: 



L\NH 



1 — Pr( System empty if no HP packet) = 1 



P^ 
P^' 
(80) 



As Pq is known (Eq. [50]), we have two equations (Eq. [79] and 
80) and two unknowns (E[S] and Pl\nh)^ so they can be 
Tound. 

For the LP packets, the waiting time can be found easily 
from the conservation law. To find the completion time, the 
same approach which was used as an approximation for small 
availability duration (Eq. [73]) can be used here because the 
exponentially distributed service time of LP packets makes 
the process for LP packets a renewal process. In other words, 
the distribution of interruptions, from the LP packets point of 
view, is identical. The only change is that here Y should be 
replaced with y<^2 _ Y\{Y < T2) because the completion 
time continues if the availability period is smaller than the 
remaining service time of the LP packet (T2). We have: 



p . 

->- at 



Ai 



aiY<T 



Ai +a + 7' 



PaiR = (1 - PaiY<T)Pr{R < Ai), 
PNaiC< = (1 - PaiY<T - PaiR) 

a + 7 



(81) 



(82) 



E[B, 



l-\iE[Xl] 

E[Rr]+E[Xl] 
^[^^"l l-\,E[Xl] ' 



{l-Pr{R<A^)), (83) 
(84) 



Ai + a + 7 

E[R]^^^E[Xl] 



(85) 



E[R] = PNa^C<{E[R<^^]) + Pa^Y<T{E[Bc^ - E[Y<^]) 
-^Pa^R{E[Bn^]^E[A<'']). 

(86) 

Similarly, E[R2'^] and consequently E[X2] and E[X2] can 
be found. 

It is worth noting that for two classes of traffic, only the 
assumption of Exponential service time for LP packets is 
required to find the results above. The service time of HP 
packets can be general. 




0.04 0.05 0.06 

LP Arrival Rate {X^ 



Fig. 8. System time of HP packets vs. LP arrival rate, for small and large 
scenarios. 



IV. Simulation Results 

In this section, simulation results for two different scenarios 
are provided: when interruption durations are distributed ex- 
ponentially and when they are constant. In both scenarios, we 
assume two classes of traffic: high priority (HP) and low prior- 
ity (LP), respectively known also as type-1 and type-2 pack- 
ets. The service time (packet length) is also assumed either 
Exponential or constant. The scenarios are called 'ExpExp', 
'ExpDet', 'DetExp' and 'DetDet' in the figures (respectively 
for T and R). Four disciplines discussed are called respectively 
'Non' (non-preemptive), 'ENo' (exceptional non-preemptive), 
'Pr' (preemptive) and 'FP' (preemption in case of failure) in 
the figures. 

The duration of availability periods is selected in a way that 
it models realistic scenarios: when it is much larger than the 
service time of a single packet (called large scenario where 
E[Y] = 75 and E[R] = 15), which represents an almost static 
or lowly dynamic CR network, and when it is much smaller 
than the service time which represents a highly dynamic 
(opportunistic) CR network |T6| (called small scenario where 
E[Y] = 1 and E[R] = 0.2). Their ratio. Pa, is selected the 
same in both scenarios. 

We assumed that HP packets model signaling and control 
packets. Thus, their arrival rate, Ai is usually lower than 
LP packets and the packets are also smaller. HP arrival rate 
is assumed constant (unless mentioned) equal to 0.03 and 
we have E[Ti] = 3 and E[T2] = 5. To model reahstic 
scenarios, the unit of time can be assumed milisecond (ms) 
which represents, for instance, 378-Byte and 625-Byte packets 
over a link with a data rate equal to 1Mbps. 



A. Exponential Recovery Time 

Exponential recovery time is a realistic assumption when 
channel selection is done randomly, by sensing a list of 
channels one by one, and the user stops after finding the first 
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Fig. 9. System time of LP packets vs. LP arrival rate, for small and large Fig. 10. System time of HP and LP packets vs. the variation of availability 
scenarios. periods' duration (large scenarios). 



channel. The remaining time of the busy period after an arrival 
in R will be still Exponential. So, Rr = R. 

As expected and can be verified in all figures from now 
on, the system time for any class of traffic using FP and 
ENo disciplines are between two other schemes, and ENo 
is generally close to non-preemptive (Non). As illustrated 
in Figures [8] and|9] when the packets length is Exponential 
(ExpExp), Eq. ( fTSJ ) can be used to find the waiting time of 
HP packets. Then, Conservation Law (CL) has been used to 
find the waiting time and consequently the completion time of 
LP packets for the FP scheme. Results are represented both for 
small and large availability periods. As expected, the system 
time of HP packets in preemptive scheme is independent of 
the arrival of LP packets, so it is constant. Moreover, we 
can see in both figures that with the same ratio of E[Y] 
and E[R], performance metrics are worse for all classes of 
traffic when availability and unavailability periods compared 
to packet length is large (large scenarios). The other point that 
can be verified is the performance difference of the exceptional 
non-preemptive (ENo) and preemptive in case of failure (FP) 
compared to other two schemes. In small scenarios, the service 
can be preempted very soon because interruptions are frequent. 
The FP discipline performance is thus close to the preemptive 
scheme. Toward large scenarios, FP performance gets farther 
from the Pr discipline and gets closer to non-preemptive 
discipline. 

The ENo discipline performance is normally very close 
to non-preemptive scheme performance in small scenarios 
because the probability of an HP arrival in the first recovery 
period is very low. In large scenarios, their performance 
difference increases. 



In Figure [T0| arrival rates are assumed constant and equal 
(Ai = A2 = 0.03) and the system time of both packet types 
are shown versus the variations of the duration of availability 
periods (large scenario). The lowest value of E[Y] is selected 
equal to E[R] = 15 because we assumed that it does not 
make sense to use Opportunistic Spectrum Access (OSA) 
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Fig. 11. System time of HP and LP packets vs. the variation of recovery 
periods' duration (large scenarios). 



when availability periods are shorter than the recovery periods. 
Figure [TT] inversely shows the system time versus the variation 
of the recovery time when the average of availability periods' 
duration is assumed equal to ^[75]. 



In Figure 12 the packet lengths are both constant (DetExp). 
Thus, a lower bound for the waiting time of LP packets in the 
new preemptive scheme has been presented. Moreover, the 
results are generated for two different values of HP arrival 
rate. To see the accuracy of completion time approximations, 
simulation results and approximations for the moments of the 
completion time of LP packets, D2, are compared in the upper 
part of Table [l| Note that the completion time of LP packets 
is independent of their arrival rate. 

B. Constant Recovery Time 

Constant recovery time indicates, for instance, the scenarios 
in which the information of channels' occupancy is provided 
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Fig. 12. System time of LP packets vs. LP arrival rate for two values of HP 
arrival rate (small scenario). 
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Fig. 13. System time of LP packets vs. LP arrival rate (ExpDet, large 
scenario). 



in advance; therefore, no random sensing is required and the 
recovery time only represents a constant time for negotiation 
and radio ahgnment. In order to compare the results with the 
previous scenario, we assume the same average. The remaining 
time after an arrival in R will be distributed hy Rr = R — 
A\{R>A). 

When real service time (packet length) is Exponential 
(ExpDet), waiting time of HP and LP packets in the FP 



discipline can be found. Figure [13] illustrates the system time 
of LP packets versus their arrival rate. The system time of 
HP packets in the same scenario is sketched in Figure [T4| In 



the last scenario where it is assumed that the packet length is 
also constant (DetDet), and the results are shown in Figures 



13] and [T4| respectively for LP and HP packets. It should be 
noted that the upper bound for HP packets in this case is loose, 
so it is not sketched to enhance the clarity of the figure. As 
expected, it can be seen in both figures that the performance. 



Fig. 14. System time of HP packets vs. LP arrival rate (DetDet and ExpDet, 
large scenario). 

TABLE I 
Moments of the LP completion time with the policy of 

PREEMPTION IN CASE OF FAILURE (FP) FOR DETEXP AND DETDET 

SCENARIOS (S: Simulation, A: Approximation). 



Seen. 


Ai 


E[X2]-S 


E[X2]-A 


ElXl]-S 


ElXi].A 


Small,DE 


0.03 


6.62 


6.72 


47.01 


49.51 


Small,DE 


0.05 


7.10 


7.32 


57.05 


62.33 


Large,DE 


0.03 


6.13 


6.12 


77.53 


74.96 


Large,DE 


0.05 


6.26 


6.25 


88.37 


85.62 


Small,DD 


0.03 


6.61 


6.72 


46.76 


49.23 


Small,DD 


0.05 


7.10 


7.31 


56.69 


61.98 


Large,DD 


0.03 


6.14 


6.12 


58.40 


56.45 


Large,DD 


0.05 


6.26 


6.25 


65.00 


63.64 



when the distribution of real service time (packet length) is 
Exponential, is worse compared to the case where the packet 
length is constant (with the same average). 

Again to see the accuracy of approximations for the com- 
pletion time of LP packets, simulation and analytical results 
are compared in the lower part of Table |l| for two different 
values of HP arrival rate. 

V. Conclusion and Future Work 

Priority queueing is a classical scheme to provide traf- 
fic differentiation in communication links. To tackle various 
priority queueing schemes for opportunistic spectrum access 
(OSA), implemented by cognitive radios, we discussed in 
this paper first a general queueing model with interruptions. 
The derived results were used to mathematically solve three 
different priority disciplines in the presence of interruptions: 
non-preemptive, exceptional non-preemptive and preemptive- 
resume, and to discuss and provide some approximations for a 
new priority scheme which is preemptive in case of failure. By 
simulation, derived results and approximations were validated 
for different scenarios. 

As discussed in the introduction, one application of these 
derived results can be to use them in an optimization and 



decision-making problem in cognitive radio networks where 
the cognitive radio node is able to change the priority queueing 
scheme on the fly, and the objective is optimizing some 
performance metrics for different classes of traffic. In such a 
decision-making model, the decision variable is the priority 
queueing discipline to be selected among these proposed 
schemes. This implies that for instance when the number of 
accumulated LP packets increases, the CR node may start 
to employ the new scheme of preemption in case of failure, 
instead of a pure preemptive. Derived results is used to find 
the immediate cost of making a decision. This idea will be 
explored in our future work. Moreover, queueing with service 
repeat after an interruption and with variable service rate will 
be also discussed in future. 

Appendix 

For two arbitrary random variables X and Y, the distribu- 
tion of two new random variables, derived from X and Y, are 
used several times throughout the paper: Z = X\{X < Y) 
and Q = X — Y\{X > Y). For simplicity of presentation, we 
discuss the statistics of these two random variables here. We 
have: 



fz{t) = 



Pr{Y>t)fx{t) (l-Fy{t))fx{t) 



Pr{Y > X) 



Pr{Y > X) 



(87) 



When Y is exponentially distributed with parameter a, we 
have: 







*/x(t) 



Pr{X <Y) = 

X[a) 
In this case, E[Z] can be given by: 

^ -d/daX{a) 
^ ^ X{a) 

The second moment can be given similarly: 

d2/da'^X{a) 



E[Z^ 



X{a) 



(88) 
(89) 

(90) 

(91) 



For the second random variable, Q, we still assume that Y 
is Exponentially distributed with parameter a. Then, by doing 
algebra (details can be found in |6, Lemma2]) or in JSJ, we 

have: 

E[X] _ _ 1 

a 



E[Q] = 



1 



and 



E[Q' 



X{a) 

E[X^]-2^ 
1 - X{a) 



(92) 



a 



2 • 



(93) 
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TABLE II 
Notations 



Notation 



Description 



Y 

R 

C 

X 

A 

X 

Index jj 

Index e 

Index a 

Index ^u 

Index u 

B 

Bb 

Bz 

T 

O 

Jc 

Js 

Z{s) 

fz{t) 

Fz{t) 

m(t) 

'rna\b\u(t) 

m'^(t) 

P* 

P 

Tr 

Pa 
Pae 

w 

D 

NQ 

PaiC\Y\R 

K 

1 = 1,2 

OL 

P 
FzO- 

7 

x(t) 



Length of availability periods (RV) 

Length of recovery (interruption) periods (RV) 

C=Y+R 

Arrival rate 

Inter-arrival time 

Completion time (RV) 

For packets entered a busy system 

For packets entered an empty system 

For packets entered an empty-available system 

For packets entered an empty-unavailable system 

For packets whose service started at the beginning of a y 

HP busy periods (RV) 

Busy period started with X^ 

Busy period started with Z + X^ 

Real service time (RV) 

Min{C,T2) (RV) 

Remaining completion time of the packet in service 

Remaining service time of the packet in service 

LST of a continuous random variable Z 

P.D.F of a random variable Z 

C.D.F of a random variable Z 

Average number of renewals until time t 

m(t) for packets of type a, b or u 

Second moment of the number of renewals until time t 

A£;[x*] 

\E\Xi,] 

Remaining transmission time in the last period (no interruption) 

Probability of system being available 



E[Y] 



^r 



<Z2 



S 
Po 



E[Y]-\-E[R] 

Probability of system being available when empty (arrival point) 

Waiting time in the queue 

Total time spent in the system 

Average number of packets waiting in the queue 

Probability of HP arrival in a cycle (C), Y or R 

Local parameter to count the number of an event 

Superindex or Subindex represents the traffic class 

Exp. Dist. parameter when Fy = 1 — e"""* 

Exp. Dist. parameter when Fr = 1 — e~^* 

Defined in Eq. Jl6| as the minimum of random variable Z and an Exponential distribution 

Exp. Dist. parameter when Ft = 1 — e~'^* 

Completion time of a packet with length t over a queue with interruption 

Number of priority classes 

Service time of the virtual packets which form the interruptions 

Remaining of R after an arrival in R 

Remaining of C after an arrival in C 

Zi\(Zi < Z2) (for two RVs) 

Initial setup time (RV) (Section II-H 

Probability of system being empty 



